Orthorhombic crystalline field splitting on orbital and magnetic orders in KCuF 3 
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Magnetic and orbital structures in KCUF3 are revisited by the cluster self-consistent field approach 
developed recently. We clearly showed that due to the inherent frustration, the ground state of 
the system with the superexchange and Jahn- Teller phonon-mediated orbital couplings is highly 
degenerate without broken symmetry; the orthorhombic crystalline field splitting arising from static 
Jahn- Teller distortion stabilizes the orbital ordering, about 42% in the x 2 — y 2 orbit and 58% in the 
3z 2 — r 2 orbit in sublattices. The magnetic moment of Cu is considerably reduced to 0.49/is, and 
the magnetic coupling strengths are highly anisotropic, Jc/Jab ~ 26. These results are in agreement 
with the experiments, implying that as an orbital selector, the crystalline field plays an essential 
role in stabilizing the ground state of KCUF3. The ls-3d resonant X-ray scattering amplitudes in 
KCUF3 with the type-a and type-d structures are also presented. 
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I. INTRODUCTION 

Three-dimensional pseudocubic perovskite KCUF3 has 
attracted extensive interest since 1970s for its unusual 
low-dimensional antiferromagnetic (AFM) and orbital 
ground state. In this compound the Cu 2+ ion has 3d 9 
configuration with the fulfilled t^g orbits and the twofold- 
degenerate e g orbits occupied by one hole, the latter leads 
to the Jahn- Teller (JT) effect. The orbital degree of free- 
dom of the holes interplays with spin and lattice degrees 
of freedom, which results in the orbital regular occupa- 
tion in real space, i.e. the orbital ordering. The orbital 
polarization of the hole is usually depicted by a pseu- 
dospin operator, f = 1/2 ^ ab c\a a bCb, here c\ creates a 
hole at orbit a, and a denotes the Pauli matrix. r 2 z = 1/2 
represents the full orbital polarization in |3z 2 — r 2 ) and 
t[- = —1/2 the orbital polarization in \x — y 2 ). The or- 
bital polarization degree represented by (tj) is called the 
orbitalization. Up to date, the orbital ordering is consid- 
ered as the essential factor in stabilizing the abnormal 
magnetic structure in KCUF3, and is observed recently 
by the resonant x-ray scattering (RXS) experiment* 1 ^. 
Though it has been known very early that the combi- 
nation of the electronic superexchange (SE) interaction 
and the JT effect is responsible for the orbital ordering 
in KCuFg^, as shown in the following, the uniqueness 
of the orbital ground state and the magnetic and the 
orbital experimental results available are not well under- 
stood theoretically 3 : 4 : 5 : 6 : 7 : 8 . 

In pseudocubic perovskite KCuF 3 the adjacent CuF 6 
octahedra in (a,b) plane are elongated along the a and 
b axis alternatively due to the JT distortion. There are 
two tetragonal crystal polytypes discriminated as type- 
a by the antiferro-distortion stacking the ab-planes and 
type-d by the ferro-distortion stacking the ab-planes&iii. 
The RXS results^ showed that in KCuF 3 the ground 
state is the G-type antiferro-orbital (AFO) order for type- 
a structure; and the C-type AFO arrangement for type-d 



structure. In these two structures the Cu spins are A- 
type AFM order below 39 K for type-a and below 22K 
for type-d; in both situations the averaged magnetic mo- 
ment of each Cu ion is about 0.49 ijlb at 4 Kii. The mag- 
netic couplings in KCUF3 are highly anisotropic, and the 
neutron scattering experimenljAiiS gave rise to J a fc/J c ~ 
-0.01, suggesting significant one-dimensional character. 

More than thirty years ago Kugel and Khomski: 3 pro- 
posed a SE model from the pure electron-electron inter- 
action, the so-called K-K model, between Cu 3d elec- 
trons to describe the role of orbital degree of freedom in 
the A-type AFM structure; by comparing the energies of 
different magnetic structures in the classical approxima- 
tion, they found the A-type AFM is stable; while for the 
orbital configuration, the G-AFO ordering is degenerate 
with the C-type AFO ordering. Further they showed^ 
that the anharmonic distortion from the JT effect low- 
ering the lattice symmetry might give rise to the orbital 
ordering. In fact the orbital part of the K-K model is in- 
herently frustrated^; Feiner et ulJ±H pointed out that 
due to quantum spin-orbital wave excitations, this frus- 
trated SE interaction leads to a spin-orbital liquid state; 
and when a tetragonal crystalline field (CF) splitting, 
E 2 , is applied to the K-K model, the ground state is 
ferro-orbital at large E z . Some other authors 7 , empha- 
sized that the phonon-mediated orbital coupling arising 
from the cooperative JT effect gives rise to the orbital 
ordering. However, Khomskii and Mostovoy 1 ^ recently 
pointed out that this effective orbital-orbital interaction 
is also inherently frustrated, similar to the orbital part 
in the electronic SE interaction. In this case the JT or- 
bital interaction combining the SE interaction favors the 
para-orbital or orbital liquid phase since the quantum 
fluctuations of the pseudospins f are still very large. 

Under the cooperative JT distortion, the low- 
temperature crystal structure of KCuF 3 is tetragonal; 
meanwhile, the local crystalline field of GvFq arising 
from the static Jahn- Teller distortion is orthorhombic. 
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The effect of the orthorhombic CF splitting on the or- 
bital order and the magnetic order of KCuF 3 was seldom 
taken into account in the past literatures. In the recent 
ab initio study within the LDA+U scheme, Binggeli and 
Altarelliii found that in the type-a JT distorted tetrag- 
onal structure, the orbital ordering is G-type AFO and 
the orbital RXS intensity agrees with the experimental 
observation; in contrast, Medvedeva et al& found in the 
absence of the JT distortion, the C-type AFO ground 
state is degenerate with G-type AFO. Therefore the re- 
alistic orthorhombic CF splitting is crucial for the stable 
orbital-ordered ground state. Nevertheless, the ab initio 
study underestimates the spin and orbital quantum fluc- 
tuations in KCUF3. It deserves to explore in detail on 
how the spin and orbital quantum fluctuations and the 
role of orthorhombic CF splitting on the orbital ordered 
ground state. 

In this paper, after deriving the orthorhombic CF 
splitting, we study the combination effect of the elec- 
tronic SE coupling, the effective orbital JT coupling and 
the CF splitting on the orbital and spin ground state. 
Utilizing the cluster self-consistent field (Cluster-SCF) 
approach 20 , we demonstrated that the orthorhombic CF 
splitting plays a key role in stabilizing the orbital or- 
dered phases of KCUF3 in type a and type d structures: 
the orbital and spin fluctuations considerably reduce the 
magnetic moment of Cu spin to 0.49/is, and the strong 
anisotropy in spin correlation functions and orbital cor- 
relation functions results in the ratio of magnetic cou- 
pling strengths, J c /J a fc ~ 26. The azimuthal dependence 
of the RXS intensity is also calculated for the Is — 3d 
excitation event. The experimental results could be con- 
sistently understood in the present theory. The rest of 
this paper is organized as follows: in Sec. II we describe 
the effective Hamiltonian and the Cluster-SCF method; 
then we present the theoretical results and discuss the 
role of CF in magnetic and orbital orderings in Sec. Ill; 
the azimuthal angle dependence of the RXS intensity is 
given in Sec. IV; and the last section is devoted to the 
remarks and summary. 

II. MODEL HAMILTONIAN AND 
CLUSTER-SCF METHOD 

According to the preceding analysis, the effective 
Hamiltonian of KCUF3 contains three parts 

H = Hse + Hjt + Hcf (1) 

The first term Hse represents the electronic SE cou- 
pling between two nearest-neighbor (N.N) e g holes of 
Cu 2+ ions derived from the generalized twofold degen- 
erate Hubbard modeU^, 

H S e = ( Jl ^ ■ *3i + J 2 J M ' % + Jslll^Si ■ s jt 

l=x,y,z 

+ J 4 /|4) (2) 



where the constants Ji, J2, J3 and J4 
are the superexchange coupling strengths: 
Ji = 8t 2 [U/(U 2 -J 2 H )-J H /(U*-J H )], 

j 2 = m 2 [1/(11! + j h ) + i/(u + j h )}, j 3 = 

32t 2 [U 1 /(U 2 - J 2 H )- J H /(U 2 - J%)], and J 4 = 
8t 2 [{U x + 2 J H )/(U 2 - J%) + J H /(U 2 - J%)] , respec- 
tively. U and U\ are the intra- and inter-orbital Coulomb 
interactions, and Jh is the Hund's coupling. Due to the 
pd hybridization between Cu 3d and F 2p orbitalsi£, we 
take the relationship U = U\ + Jh- The parameters 
U=7.5 eV and J#=0.9 eV are adopted from the con- 
strained LDA computation for KCuFa 4 . The hopping 
integral along the c direction is the largest, t3 Z 2_ r 2 3z 2_ r 2 
= At, we take t= 0.12 eV, thus the energy scale of the 
superexchange coupling is J = I6t 2 /U = 30.7 meV. 
Si denotes the spin at site i, while the operator 
l\ = cos (27TTOi/3) rf — sin (27rm;/3) rf , the index I de- 
notes the direction of a bond, I = x, y or z, corresponding 
to the crystal axes, a, b and c; (ij)i connects site i and 
its nearest-neighbor site j along the I direction, and 
(m x ,my,m z )=(l,2,3). The orbital pseudospin operators 
r? and rf are the components of f. 

The second term in Eq.(l) represents the phonon- 
mcdiatcd orbital coupling from the cooperative JT effect. 
Through eliminating the phonon operator, one can get an 
effective orbital-orbital interaction Hj j V 2 i 19 ; 

H JT = gjT ll A ( 3 ) 

l=x,y,z 

It has the same form as the orbital part of the last term in 
Eq.(2), both of them contribute to the orbital frustration 
and quantum fluctuations. The Jahn- Teller distortion 
energy Ejt ~ 130 meV—, gjr is the same order as Ejt 
in magnitude. 

The static JT effect distorts the symmetry of CuFg oc- 
tahedra to orthorhombic, in which there are three differ- 
ent Cu-F bonds, the medium length Cu-F bond is along 
the z-axis, the long and short Cu-F bonds alternate along 
the x and the y axes. All of the F-Cu-F angles are 90° 
or 180°, and no rotation is found* 10 . Employing the point 
charge model to calculate the CF splitting, we obtain the 
orthorhombic CF splitting of the holes, 

H C F = J2( V ^i +Vi*T?) (4) 

i 

The intermediate compression of the Cu-F bond along c- 
axis implies the sign of Vi Z is negative for the hole, Vi Z < 
0, since the compressed octahedra lifts the |3z 2 — r 2 ) 
orbit, thus the |3z 2 — r 2 ) orbit is in favor of hole occu- 
pation; and the component V% x mixes the two e g orbitals 
at site i. The former favors of the orbital ordering, while 
the latter tends to destroy the orbital ordering. Using 
the crystal structure data in Ref.10 we estimate the ratio 
\V X /V Z \ ~ 2, and V z is the same order in magnitude as 
the superexchange interaction. 

It is a huge challenge to treat the spin-orbital correla- 
tions and fluctuations and to find the ground state of such 
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strongly correlated systems with high accuracy. To study 
the groundstate properties, we apply the Cluster-SCF 
approach^!! developed recently to deal with the compli- 
cated spin-orbital Hamiltonian (1). This approach com- 
bines the exact diagonalization for a central cluster and 
the self-consistent field for surrounding atoms. The main 
idea of the approach is described as follows: consider a 
proper cluster, usually the unit cell of the compounds, 
in which the 3d electrons interact via Eq.(l). First, we 
substitute the spin coupling s, ■ Sj in the cluster Hamil- 
tonian, H, with the spin correlation functions (s, • Sj), 
here 

H = E [ F tt> 3i ) $ ' g ii + ( J 4 + 9jt)I\ 4 ] 

l=x,y,z 

+ ^2(V iz Tf + V ilc 7f), (5) 
i 

with 

F (Ti, fj) = Jl + y {I\ + 7j, ) + J 3 ^4 ■ 

And then diagonalize the orbital part of the cluster 
Hamiltonian in the presence of the orbital SCF 21 , hence 
obtain the orbitalization (r) and the orbital correlation 
functions (7$ • fj). Second, substitute the orbital oper- 
ator, r, and their coupling fi ■ fj with the orbitaliza- 
tion and the orbital correlation functions obtained, and 
diagonalize the spin part of the cluster Hamiltonian in 
the presence of the spin SCF, thus obtain a new set of 
spin correlation functions. Repeat the above steps it- 
eratively until the groundstate energy, the spin and the 
orbital correlation functions self-consistcntly converge to 
the accuracies. The advantage of the present approach 
superior to the traditional mean-field method is that the 
short-range spin and orbital correlations and quantum 
fluctuations are taken into account. In comparison with 
the traditional mean-field method, one can obtain more 
better results for some simple models with large quantum 
fluctuations, such as the Heisenberg AFM modelSi. 

III. RESULTS AND DISCUSSIONS 

In this section we first investigate the ground state 
of KCUF3 under the electronic SE interaction and the 
JT phonon-mediated orbital interaction, then explore the 
role of orthorhombic CF splitting in the ground state. 

A. Superexchange and JT Orbital Interactions 

For the electronic SE coupling in Eq.(2), the mean- 
field results^ suggested the staggered orbital order with 
x 2 -z 2 / y 2 -z 2 orbits. However, utilizing the cluster-SCF 
approach, we find the spin-orbital ground state is com- 
posed of numerous degenerate states: such as the Neel 



AFM order with ferro-orbital structure, in which the hole 
occupies one of the three orbits among the \3x 2 ~r 2 ), the 
|3y 2 - r 2 ) and the |3z 2 - r 2 ) orbits; or the Neel AFM or- 
der with the alternating plaquette valence-bond order, 
in which one plane, for example, the xy-plane, is occu- 
pied by the |3x 2 — r 2 ) orbit and the nearest neighbor 
ones is taken up with the \iy 2 — r 2 ) orbit, etc. The de- 
generate ground state contains the alternating plaque- 
tte valence-bond component, in agreement with Feiner 
et al.'s result^ by considering Gaussian quantum fluc- 
tuations. This also confirms the validity and the effi- 
ciency of our Cluster-SCF approach. In this case, the 
A-type AFM structure observed in realistic KCUF3 is 
not the candidate of the ground state. On the other 
hand, the nearest-neighbor orbital correlation functions 
are significantly different from zero. The short-range or- 
bital correlations are strong, indicating that the system 
with only the electronic SE interaction is a spin-orbital 
liquid phase. 

In such a spin-orbital correlated system, the spin align- 
ment strongly depends on the orbital configuration. In 
fact one could easily find that the system with the SE in- 
teraction is a spin-orbital liquid phase, rather than a spin- 
orbital ordered state. Due to the spin rotation symmetry, 
the spin interaction in the SE coupling is Heisenberg-like. 
After averaging over orbital freedom degree, the spin ex- 
change coupling strength along the Z-axis in Eq.(5) reads 

= Ji + y «4> + <4» + J 3</i/j>, 

which is isotropic for the x-, the y- and the z-axes, such 
an interaction does not lead to the anisotropic A-type 
AFM structure in KCUF3, unless the orbital symmetry 
is broken. In the orbital part of the SE interaction in 
Eq.(2), after averaging over the spin coupling, the orbital 
part along the Z-direction reads, 

h l = y (Si ■ s h ) (I\ + ij,) + (Ms, ■ 4) + J 4 ) III], , (6) 

The first term is similar to a magnetic field, called the or- 
bital field. This orbital field is frustrated for the orbital 
occupation, which favors either \iz 2 — r 2 ) or \x 2 — y 2 ) 
orbits along the z-axis; while along the x direction, it 
favors either |3x 2 — r 2 ) or \y 2 — z 2 ) orbits. Unfortunately 
the orbital interaction part, /|/j £ , in Eq.(6) is also in- 
herently frustrated: no matter what the coefficient of the 
second term in Eq.(6) is, negative or positive, this orbital 
correlation favors a frustrated ground state, as pointed 
by Khomskii et al^ On the one hand, if the coefficient 
of jjjj is negative, the orbital exchange coupling favors 
different ferro-orbital occupations along different direc- 
tions, in accordance with the frustration of the orbital 
field mentioned above; on the other hand, if the coeffi- 
cient is positive, the orbital polarizations in each bond 
arising from the orbital field are opposite to that arising 
from the orbital coupling, which enhances the frustration 
effect. 
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One notices that if the spin correlation functions are 
strongly anisotropic, the frustration effect is greatly sup- 
pressed. For example, as the spin correlation along the 
z-axis is so strong that (s^ • Sj z ) rs —3/4, while the spin 
correlation functions along the x, y-axes almost vanish, 
(si ■ Sj^ ) rs 0, then the orbital part of the SE interaction 
becomes 

h = ./-.^ /;/;. , + Iflf +y - 0.95I?Il +z - 1.78/?] (7) 

i 

the x, y-components of the orbital field approach to zero, 
leaving a large z-component in Eq.(7). Obviously the 
strong uniaxial orbital field suppresses the frustration 
of the orbital exchange coupling, and singles out the 
\iz 2 — r 2 ) orbit in each site, forming the ferro-orbital or- 
der; our numerical calculation confirms this result. Also 
it could be shown that strongly anisotropic plaquette- 
valence-bond correlation of the spins favors the alter- 
nating plaquette-valence-bond orbital order. As we will 
shown later, the highly anisotropic magnetic correlation 
is the consequence of the orbital ordering in KCuF 3 . 

Neglecting the spin and orbital quantum fluctuations, 
Kugel and Khomskii? found the mean-field solution of 
the SE interaction is the A- type AFM and G-type/C- 
type AFO order. In fact, their classical approximation 
to the spins as the A-type AFM order introduced an or- 
bital field li)- This orbital field lifts most degener- 
ate states of the orbital ground state, leaving the G-type 
and the C-type AFO configurations as the candidates; 
the fourfold rotation symmetry of the spin structure also 
confines the possible orbital ground state as the G-type or 
C-type AFO structure, in agreement with Goodenough- 
Kanamori empirical rule 2 ^. Nevertheless, the orbital field 
do not remove the degeneracy of the G-type and the C- 
type AFO structures, as shown in Ref.3. However, as 
found in the experiments^, the average magnetic mo- 
ment of each Cu spin is about 0.49 /is, only a half of 
the classical expectation, suggesting that spin quantum 
fluctuations in KCUF3 is very strong, and the classical 
approximation to the spins is not appropriate. 

The phonon-mediated orbital-orbital coupling with 
discrete cubic symmetry, Hjt in Eq.(3), is also strongly 
frustrated. It has the same impact on the orbital ground 
state as the second term in Eq.(6) does. Early treatment^ 
to Hjt by the mean-field approach is questionable since 
the presumed orbital order was equivalent to introduce 
an artificial orbital field to break the orbital symmetry, 
however the orbital field did not really exist. We notice 
that Hjr can be absorbed in the second term in Eq.(6), 
so the combination of the JT orbital coupling Hjx and 
the SE interaction, Use, does not break the symmetry 
of spin and orbital in cubic crystal structure. Therefore 
in the absence of the crystalline field splitting from the 
static JT distortion, the ground state of the system is an 
orbital liquid or para-orbital phase. Accordingly, the or- 
bital symmetry in KCUF3 should be broken by the static 
JT distortion, i.e. the orthorhombic CF splitting, Hcf- 
In what follows we explore the role of orthorhombic CF 



splitting in the ground state in KCUF3. 



B. Role of Orthorhombic Crystalline Field 

In perovskite KCUF3 there exist two kind cooperative 
JT distortions at low temperatur e 9 ! 10 . With respect to 
these two different small lattice distortions, KCUF3 ex- 
hibits two slightly different crystalline phases, the type-d 
and the type-a structures. The orbital ground state in 
the type-a structure differs from that in the type-d struc- 
ture due to different orthorhombic CF splittings of the 
CuFg octahedra. 

We first study the orbital components of single Cu 
atom under the orthorhombic CF. For |T4/T4| = 2, the 
orbital wavefunctions always consist of two orbital pat- 
terns, a low energy pattern \a) — 0.851|3z 2 — r 2 ) ± 
0.526|3x 2 — y 2 ), in which the hole completely occupies 
this pattern at large | V z \ , and a high energy pattern 
1 6) = 0.526|x 2 - y 2 ) =F 0.851|3z 2 - r 2 ); here '±' refer 
to the two sublattices of the antiferro-distortion in the 
two crystalline phases. Approximately, \a) ~ \y 2 — z 2 ) 
and \b) w \ix 2 - r 2 ) for '+'; and \a) w \x 2 - z 2 ) and 
1 6) « \3y 2 — r 2 ) for Such combinations are also con- 
sistent with the orthorhombic distortions of the CuFg 
octahedra: if the Cu-F bond is elongated along the x 
axis, the CF singles out the \y 2 — z 2 ) orbit, correspond- 
ing to '+'; on the other hand, if Cu-F bond is elongated 
along y direction, the energy of \x 2 — z 2 ) orbit is lower, 
corresponding to 

In the lattice case, the orthorhombic CF splitting com- 
petes with the SE interaction and the JT orbital coupling, 
the orbital occupation of Cu 3d holes depends not only 
on the CF splitting ratio 114/141, but also on the mag- 
nitude of V z . When \V Z \ is very small, due to the orbital 
frustration and large quantum fluctuations, the orbital 
symmetry is not broken, and the ground state of the sys- 
tem is still an orbital liquid or para-orbital phase. The 
critical value of \V Z \ breaking the orbital symmetry relies 
on the JT orbital coupling. For large |14|, the transverse 
CF splitting V x alternating in the xy-plane in the type-d 
structure gives rise to C-type AFO configuration; as a 
contrast in the type-a structure, the staggered transverse 
CF splitting in the x, y and z directions gives rise 
to G-type AFO configuration. At IT4H0.5J, the sublat- 
tice orbitalization and the orbital correlation functions 
listed in Table I definitely show the G-type AFO corre- 
lation in the type-a structure and the C-type AFO corre- 
lation in the type-d structure. Some magnetic properties 
in both structures are also collected in Table I. Therefore 
the presence of the orthorhombic CF breaks the discrete 
orbital symmetry, suppresses the orbital frustration and 
quantum fluctuation, and establishes the long-range or- 
bital order in KCUF3. 

As soon as the orthorhombic CF singles out the or- 
bital structure, it also stabilizes the magnetic structure 
simultaneously, hence the spin-orbital ground state. Un- 
der the full Hamiltonian, our numerical results show that 
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the magnetic ground states are the A-type AFM order 
both for the type-a and for the type-d crystalline phases. 
And the spin correlations are strongly anisotropic, (s*j • 
Sj z )/(si ■ Sj x ) rs 10. Furthermore, we obtain the spin 
coupling strengths J z and i x ,y, which are 16.6 meV and 
0.64 meV, respectively, giving rise to \ J Z /J x ,y\ about 26 
for V Z =-0.5J, as shown in Fig.l. These results are in 




FIG. 1: Dependence of the ratio of N.N spin coupling 
strengths J x ,y/Jz on orthorhombic CF splitting |14|/J. Here 
J is the energy scale of the SE coupling. 

agreement with the experimental data in K(.'uF : < l 2 1 1 . 
Such strong anisotropy in spin correlations and magnetic 
couplings is attributed to the anisotropic regular distri- 
bution of the orbital wavefunctions in real space, i.e., the 
orbital ordering, as we will show later. 

Next we focus on the magnetic moment of each Cu 
spin. We find that the effective magnetic moment of 
each Cu spin is considerable reduced from 1 \ib to 0.496 
[Lb , the averaged spin of each Cu ion decreases to 0.248, 
about a half of 1/2, consisting with the neutron scatter- 
ing experimental data very wellii. The theoretical and 
experimental results are listed in Table 1. Obviously such 
great reduction of the magnetic moment arises from the 
spin-orbital quantum fluctuations: one possibility is from 
the low-dimensional AFM spin-wave excitation; the an- 
other is from the spin-and-orbital wave excitation in the 
spin-orbital systemia, which causes more spin flipping 
via the spin-orbital interaction. From Table 1, one finds 
the discrepancy between the present ratio of | J z / J x ,y\ (~ 
26) and Satija's experimental fitting data {\J Z /J x , y \ ~ 
100)2^, we attribute this discrepancy to the frozen of the 
orbital excitations in his fitting to the experimental data. 
According to the spin exchange couplings shown in Ta- 
ble I, we find that in the mean-field approximation, the 
theoretical Neel temperature of KCUF3 is about, 37 K, 
in agreement with the experimental data 39 K in type-a 
structure, confirming that our choice to the theoretical 
parameters t, U, in, gjT and V z is appropriate. 

In fact the influence of the orthorhombic CF splitting 
and the JT orbital coupling on the magnetic moment 
is not monotonously. As shown in Fig. 2, the phonon- 
mediated JT orbital coupling and the orthorhombic CF 
distortion play distinct roles in the magnetic moments 
through affecting the orbital quantum fluctuations and 



TABLE I: Calculated and experimental averaged spin, mag- 
netic couplings, nearest-neighbor spin correlations, orbital- 
ization and orbital correlations, 'a' and 'd' refer type-a and 
type-d phases. Jx,y,z are in units of meV. The theoretical 
parameters are gjr = 7J, and \V Z \ = 0.5J. 
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the orbital ordering. At small splitting \V Z \ and in the 
absence of the JT coupling, gjT = 0, the orbital field 
—If in the A-type AFM structure is much stronger than 
that from the CF splitting, resulting in large orbital po- 
larization with dominant |3z 2 — r 2 ) orbit; with the in- 
crease of the CF splitting, the transverse CF term V x 



mixes the \3z — 



orbit with the 



y ) orbit, lead- 



ing to the descent of the sublattice orbitalization, as seen 
in Fig. 2a. Meanwhile the decrease of the orbitalization 
weakens the low-dimensionality of the spin correlations, 
the magnetic moment gradually lifts with increasing CF 
splitting, which can be seen in Fig. 2b. At sufficient large 
\V Z \, the orbital occupation is full polarized at the \a) 
orbit, and the magnetic moment saturates to 0.61 [Ib 
per site. In this situation, a small CF splitting favoring 
of the A-type AFM structure produces strong local or- 
bital field and orbital polarization, suppresses the orbital 
frustration from the superexchange coupling. With the 
further increase of the CF splitting, the transverse com- 
ponent of the CF splitting mixes the orbitals |3z 2 — r 2 ) 
and \x 2 — y 2 ), the orbital polarization declines and the 
anisotropy of orbital and spin correlations become weak, 
hence the ratio of the magnetic couplings J z /J XjV de- 
crease with the increase for very large CF splitting, as 
seen in Fig.l. 

When the JT orbital coupling is taken into account, the 
dependences of the spin coupling, the sublattice orbital- 
ization and the magnetic moment on \V Z \ are different 
for gjT/J=l and for gjr/J =7- The different behaviors 
arise from the distinct effects of the CF splitting on the 
local orbital field and the frustration term. At gjy/J=l, 
\J X /Jz\ is larger than that at gjr=0, implying the low- 
dimensional characters of the spin correlations and the 
spin fluctuations become weak. This leads to large mag- 
netic moment and small anisotropy, as we find in Fig.l 
and Fig. 2b. The weakness of the anisotropy of the mag- 
netic couplings at gjT/J=l arises from the frustration en- 
hancement comparison with that at gj*r=0. Our numer- 
ical results show that at g,/T/J=l, the further increase 
in the CF splitting is almost balanced by the orbital field 
and the frustrated orbital-orbital couplings, hence the 
spin and the orbital correlations, the sublattice orbital- 
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FIG. 2: Dependence of averaged spin (a) and sublattice or- 
bitalization (b) on orthorhombic CF splitting in different JT 
orbital coupling. 



IV. RESONANT X-RAY SCATTERING FOR 
ORBITAL ORDER 

The orbital ordering in KCUF3 can be manifested in 
the RXS peaks utilizing the sensitivity of x-ray scattering 
to the anisotropic density of orbital ordered electrons^. 
The anisotropy of spatial electronic clouds gives rise to 
the anomalous tensor component in atomic scattering 
factor, and the interference of these atomic scattering 
amplitudes in the presence of long-range orbital order 
leads to the orbital superlattice reflection at the struc- 
tural forbidden position. Obviously, it is more directly 
to identify the orbital ordering using the spectral line 
shapes of the quadrupole Is — 3c? scattering, contrary to 
the complicated spectra of the Is — 4p dipole scattering 
which is often used in present experiments^, although 
the signal enhancement of the quadrupole scattering is 
weaker than the dipole scattering2i2£i2ii. In the follow- 
ing we present the azimuthal angle dependence of the 
Is — 3d RXS intensity to directly demonstrate the char- 
acter of the orbital order in KC11F3. 

For the C-type AFO ordered ground state with type-d 
structure, the orbital ordering peaks reflect at (h, k, I) — 
(odd, odd, even) , which are forbidden for the structural 
and magnetic reflections. The sublattice orbital wave- 
functions consist of two different components: \ipi) = 
ai|3z 2 - r 2 ) + a 2 \x 2 - y 2 ) and \ip 2 ) = ai\iz 2 - r 2 ) - 
ct 2 \x 2 — y 2 ), here the coefficients 011,2 are the functions of 
the interaction parameters. Then the orbital structural 
factor is read as: 



ization and sublattice magnetization and the magnetic 
coupling almost do not change with the increase of the 
CF splitting, as shown in Fig.l and Fig. 2. 

In contrast, the large JT orbital coupling at gjx = 7 J 
leads to strong orbital fluctuation, hence to small sublat- 
tice orbitalization and magnetization in Fig. 2. In this sit- 
uation, one would expect more weak anisotropy or large 
ratio \Jx/Jz\', however as shown in Fig.l, one curiously 
finds that the anisotropy is the most strong, in compar- 
ison with gjT = and gjT = J- We find this strong 
anisotropy may come accidentally in small orthorhom- 
bic CF splitting: in the spin coupling along the x-axis, 
the sign of the constant part Ji is contrary to that from 
the other terms, which leads to 3 X be very small, hence 
small ratio \J X /J Z \ and large anisotropy. Strong orbital 
fluctuations also excite spin flipping excitation via the 
spin-orbital coupling. Thus the magnetic moment crit- 
ically decreases to 0.49 fiB at \V Z \/J = 0.5. Further 
increase of the orthorhombic CF splitting greatly sup- 
presses the orbital and spin fluctuations, and the sublat- 
tice magnetic moment and the orbitalization rise a lot. 
With the further increase of the orthorhombic CF split- 
ting, the anisotropy becomes more weak, and the sub- 
lattice magnetization (S) and orbitalization (T z ) further 
increase gradually. 



Fhki = f (I\ r 2 ,ds, c, u) y/n tk n t > k iaia 2 {<izkz 

e'yk'y) + (e x k x - e y k y ) e z k' z } (8) 



(^x^x £ y^y^j 



where the function / (r, r 2 ^ s , c, ui) is the coefficient de- 
pending on the lifetime of the intermediate states, T, the 
radial matrix element r 2 ^ s , the velocity of photon c and 
the incoming photon frequency w; n e (e'),k(k') is the den- 
sity of the incoming (outgoing) beam of photons with po- 
larization e(e') and wavevector k(k'). The azimuthal an- 
gle dependence of the RXS intensity is shown in Fig. 3 at 
(1, 1,0) reflection for unrotated (era 1 ) and rotated (oV) 
channels for the perfect a polarized incoming beam. As a 
comparison, we also present the azimuthal angle depen- 
dence of the RXS intensity in Fig. 4 for the G-type AFO 
order with type-a structure at (3, 3, 1) reflection. 

Due to the difference of the orbital superlattice for the 
C-type and the G-type AFO orders, the RXS peaks ap- 
pear at different orbital reflections with (odd, odd, even) 
for the C-type orbital order in the type-d structure, 
and with (odd, odd, odd) for the G-type orbital order in 
the type-a structure. These sublattice reflections distin- 
guish these orbital orders in different crystalline phases 
of KCUF3. The azimuthal angle dependence of the RXS 
intensities exhibits different periods in the type-d and the 
type-a structures, as seen in Fig. 3 and Fig. 4. The periods 
are it and 27r for the a— ~k' channel in these two structures, 
respectively; and are tt/2 and tt for the a — a' channel 
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FIG. 3: The azimuthal angle dependence of the orbital (1,1,0) 
reflection intensity for unrotated (<xer') and rotated (<X7r') 
channels in C-type AFO. 




Azimuthal angle (degrees) 

FIG. 4: Azimuthal angle dependence of the orbital (3,3,1) re- 
flection intensity for unrotated (era') and rotated (an') chan- 
nels in G-type AFO. 

respectively. Distinct line shapes of the RXS scattering 
intensities in the a— a' channel also easily identify the two 
orbital AFO orders. The different periods are attributed 
to the different orbital structures along the z-axis. One 
may expect that the 4p orbits of Cu hybridizing with 3d 
orbits will modify the present ls-3d scattering line shape, 
but not the period. 

V. REMARKS AND SUMMARY 

We notice that in perovskite KCuF 3 , the nonresonant 
magnetic scattering experiment showed that the orbital 
angular momentum L contribute finite value the total 
magnetic moment, and L/S ~ 0.29. Obviously the finite 
orbital moment does not come from the two e g orbits, 



since the expectation of the orbital angular momentum 
L in the orbital basis wavefunctions \3z 2 — r 2 )/\x 2 — y 2 ) 
and any of their combinations is zero. One possibility of 
such considerable residual orbital moment is attributed 
to the reduced symmetry of the t2 g orbits in KCuFgSl or 
a small fraction of the ti g orbit mixing with the e g orbit; 
another possibility is from the hybridization of the 4p 
orbits with the e g orbits. In both situations the weak LS 
coupling in KCUF3 will not change the spin alignment 
considerably. 

In the numerous literatures on KCUF3, some 
authors^i^i^ took the tetragonal CF splitting into ac- 
count to stabilize the A-type AFM order, but the de- 
generacy of the G-type and the C-type AFO orders is 
not lifted, so the ground state under the tetragonal CF 
is still indefinite. Only in the present theory the full 
consideration of the orthorhombic CF splitting, V z and 
V x , together with the JT orbital coupling and the SE 
coupling, can we determine the ground state exclusively, 
and consistently interpret the experimental data. Fur- 
thermore, considering many other spin-orbital-lattice in- 
teracting compounds, such as mang aniteo 2 ?! 2 ?!??, vana- 
dium oxides22i2i, and titanium oxides22i2&Mi2£, in which 
the CF splitting arising from the lattice distortion ex- 
tensively exists, one may find such a fact that the low- 
symmetric CF splittings play crucial roles in singling out 
many degenerate candidates as the sole orbital ordered 
ground state. Thus a conclusion arrives that the highly 
degenerate ground state in the correlated electronic sys- 
tem with pure spin- orbital interactions usually stabilizes 
through distorting to a lower symmetry phase, which is a 
natural generalization of the Jahn- Teller effect in strongly 
correlated systems. Detail results will be presented in 
further study. 

In summary, we have performed a systematic study on 
the roles of the electronic SE interaction, the JT orbital 
coupling and the orthorhombic CF splitting in the orbital 
ordering and the magnetic properties in KCUF3. The SE 
and effective JT orbital coupling lead to a orbital liquid 
state due to the inherent frustration and orbital quan- 
tum fluctuations. The orthorhombic CF lowers the or- 
bital symmetry, and stabilizes the orbital ordering as the 
observed in experiment. The orbital ordering results in 
the strong magnetic anisotropy. Strong spin fluctuation 
and the orbital frustration considerably reduce magnetic 
moment of Cu spins. 



Acknowledgments 

Authors thanks N. Binggeli for providing the electronic 
structure data of KCuF 3 . Great appreciate is devoted to 
G. Sandro, M. Fabrizio and M. Altarelli in developing the 
Cluster-SCF method. Supports from the NSF of China 
and the BaiRen project from the Chinese Academy of 
Sciences (CAS) are appreciated. Part of numerical cal- 
culation was performed in CCS, HFC AS. 



8 



1 R. Caciuffo, L. Paolasini, A. Sollier, P. Ghigna, E. Pavarini, 
J. Van de. Brink, and M. Altarelli, Phys. Rev. B 65, 174425 
(2002). 

2 L. Paolasini, R. Caciuffo, A. Sollier, P. Ghigna, and M. 
Altarelli, Phys. Rev. Lett. 88, 106403 (2002). 

3 K. I. Kugel and D. I. Khomskii, Sov. Phys. JETP 37, 725 
(1973). 

4 A. I. Liechtenstein, V. I. Anisimov and J. Zaanen, Phys. 
Rev. B 52, R5467 (1995). 

5 J. E. Medvedeva, M. A. Korotin, V.I. Anisimov, and A. J. 
Freeman, Phys. Rev. B 65, 172413 (2002). 

6 J. Kanamori, J. Appl. Phys. 31, 14S. (1960). 

7 B. Halperin and R. Englman, Phys. Rev. B 3, 1698 (1971). 

8 M. Kataoka, J. Phys. Soc. Japan 70, 2353 (2001). 

9 A. Okazaki and Y. Suemune, J. Phys. Soc. Japan 16, 176 
(1961). 

10 R. H. Buttner, E. N. Maslen and N. Spadaccini, Acta 
Cryst. B 46, 131 (1990) (2001). 

11 M. T. Hutchings, E. J. Samuelsen, G. Shirane, and K. 
Yamada, Phys. Rev. 188, 919 (1969). 

12 D. I. Khomskii and M. V. Mostovoy, J. Phys. A: Math. 
Gen. 36, 9197 (2003). 

13 L. F. Feiner, A. M. Oles, and J. Zaanen, Phys. Rev. Lett. 
78, 2799 (1997). 

14 A. M. Oles, L. F. Feiner, and J. Zaanen, Phys. Rev. B. 61, 
6257 (2000). 

15 G. Khaliullin, and V. Oudovenko, Phys. Rev. B. 56, 
R14243 (1997). 

16 L. F. Feiner, A. M. Oles, and J. Zaanen, J. Phys.: Condens. 
Matter, 10, L555 (1998); M. Imada, A. Fujimori, and Y. 
Tokura, Rev. Mod. Phys. 70, 1039 (1998). 

17 N. Binggeli, and M. Altarelli, Phys. Rev. B 70, 085517 
(2004). 

18 C. Castellani, C. R. Natoli, and J. Ranninger, Phys. Rev. 
B 18, 4945 (1978). 



19 K. I. Kugel and D. I. Khomskii, Sov. Phys. Usp. 25, 231 
(1982). 

20 Liang- Jian Zou, M. Fabrizio and M. Altarelli, Preprint. 

21 Dong-Meng Chen and Liang- Jian Zou, |cond-m at/0501155 
(2005). 

22 J. B. Goodenough, Magnetism and Chemical Bond (Inter- 
science, New York, 1963) 

23 S. K. Satija, J. D. Axe, G. Shirane, H. Yoshizawa, and K. 
Hirakawa, Phys. Rev. B 21, 2001 (1980). 

24 Y. Murakami, H. Kawada, H. Kawata, M. Tanaka, T. 
Arima, Y. Moritomo and Y. Tokura, Phys. Rev. Lett. 80, 
1932 (1998). 

25 M. Fabrizio, M. Altarelli and M. Benlatlo, Phys. Rev. Lett. 
80, 3400 (1998); Phys. Rev. Lett. 81, E4030 (1998). 

26 S. W. Lovesey, K. S. Knight and E. Balcar, Phys. Rev. B 
64, 054405 (2001). 

27 M. Hidaka, T. Eguchi, and I. Yamada, J. Phys. Soc. Jpn. 
67, 2488 (1998). 

28 R. Shiina, T. Nishitani, and H. Shiba, J. Phys. Soc. Jpn. 
66, 3159 (1997). 

29 T. Mizokawa, D. I. Khomskii, and G. A. Sawatzky, Phys. 
Rev. B 60, 7309 (1999) 

30 R. Y. Gu and C. S. Ting, Phys. Rev. B 65, 214426 (2002). 

31 A. Tanaka, J. Phys. Soc. Jpn. 71, 1091 (2002). 

32 M. Cwik, T. Lorenz, J. Baier, R. Muller, G. Andre, 
F. Bouree, F. Lichtenberg, A. Freimuth, R. Schmitz, 
E. Muller-Hartmann, and M. Braden, Phys. Rev. B 68, 
060401 (2003). 

33 R. Schmitz, O. Entin-Wohlman, A. Ahaony, A. B. Harris, 
and E. Muller-Hartmann, IronTm at /0407524 (2004). 

34 M. W. Haverkort et al, Phys. Rev. Lett. 94, 056401 (2005). 

35 S. V. Streltsov, A. S. Mylnikova, A. O. Shorikov, Z. 
V. Pchelkina, D. I. Khomskii, and V. I. Anisimov, 
cond-mat/0504281 (2005). 



